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introduce the phenomenon of color superconductivity and discuss the use of weak coupling methods. 
We study the phase structure as a function of the number of flavors and their masses. We also 
introduce effective theories that describe low energy excitations at high baryon density. Finally, we 
study the possibility of kaon condensation at very large baryon density. 

Keywords. QCD, dense Matter 
PACS Nos 12.38.Aw, 12.38.Mh 

1. Color Superconductivity 

In the interior of compact stars matter is compressed to densities several times larger than 
the density of ordinary matter. Unlike the situation in relativistic heavy ion collisions, 
these conditions are maintained for essentially infinite periods of time and the material is 
quite cold. At low density quarks are confined, chiral symmetry is broken, and baryonic 
matter is described in terms of neutrons and protons as well as their excitations. At very 
large density, on the other hand, we expect that baryonic matter is described more effec- 
tively in terms of quarks rather than hadrons. As we shall see, these quarks can form new 
condensates and the phase structure of dense quark matter is quite rich. 

At very high density the natural degrees of freedom are quark excitations and holes in 
the vicinity of the Fermi surface. Since the Fermi momentum is large, asymptotic free- 
dom implies that the interaction between quasi-particles is weak. In QCD, because of the 
presence of unscreened long range gauge forces, this is not quite true. Nevertheless, we 
believe that this fact does not essentially modify the argument. We know from the theory 
of superconductivity the Fermi surface is unstable in the presence of even an arbitrarily 
weak attractive interaction. At very large density, the attraction is provided by one-gluon 
exchange between quarks in a color anti-symmetric 3 state. High density quark matter is 
therefore expected to behave as a color superconductor [1-4]. 

Color superconductivity is described by a pair condensate of the form 



Here, C is the charge conjugation matrix, and T e>,\c,tf are Dirac, color, and flavor 
matrices. Except in the case of only two colors, the order parameter cannot be a color sin- 
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glet. Color superconductivity is therefore characterized by the breakdown of color gauge 
invariance. As usual, this statement has to be interpreted with care because local gauge 
invariance cannot really be broken. Nevertheless, we can study gauge invariant conse- 
quences of a quark pair condensate, in particular the formation of a gap in the excitation 
spectrum. 

In addition to that, color superconductivity can lead to the breakdown of global sym- 
metries. We shall see that in some cases there is a gauge invariant order parameter for the 
U(l) of baryon number. This corresponds to true superfluidity and the appearance of a 
massless phonon. We shall also find that for Nf > 2 color superconductivity leads to chi- 
ral symmetry breaking. Finally, if the effects of finite quark masses are taken into account 
we find additional forms of long-range order. 

2. Phase Structure in Weak Coupling 

2.1 QCD with two flavors 

In this section we shall discuss how to use weak coupling methods in order to explore 
the phases of dense quark matter. We begin with what is usually considered to be the 
simplest case, quark matter with two degenerate flavors, up and down. Renormalization 
group arguments suggest [5,6], and explicit calculations show [7,8], that whenever possible 
quark pairs condense in an s-wave. This means that the spin wave function of the pair is 
anti-symmetric. Since the color wave function is also anti-symmetric, the Pauli principle 
requires the flavor wave function to be anti-symmetric too. This essentially determines the 
structure of the order parameter [9,10] 



This order parameter breaks the color SU (3) — > SU (2) and leads to a gap for up and down 
quarks with two out of the three colors. Chiral and isospin symmetry remain unbroken. 

We can calculate the magnitude of the gap and the condensation energy using weak 
coupling methods. In weak coupling the gap is determined by ladder diagrams with the 
one gluon exchange interaction. These diagrams can be summed using the gap equation 



Here, A(po) is the frequency dependent gap, g is the QCD coupling constant and G and 
F are the self energies of magnetic and electric gluons. This gap equation is very similar 
to the BCS gap equations that describe nuclear superfluids. The main difference is that 
because the gluon is massless, the gap equation contains a collinear cos 9 <~ 1 divergence. 
In a dense medium the collinear divergence is regularized by the gluon self energy. For 
q — > and to leading order in perturbation theory we have 
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F = 2m 2 , G=|m 2 ^, (4) 

2 |g| 

with m 2 = N jg 2 n 2 / {Air 2 ). In the electric part, m 2 D = 2m 2 is the familiar Debye screen- 
ing mass. In the magnetic part, there is no screening of static modes, but non-static modes 
are modes are dynamically screened due to Landau damping. 
We can now perform the angular integral and find 

am =4-j*to ^ ( r^s) /2 A{q :l & 

18tt 2 7 \\Po-qo\J V?o+ A (?o) 2 

with b = 256n 4 (2/Nf) 5 / 2 g~ 5 . We can now see why it was important to keep the fre- 
quency dependence of the gap. Because the collinear divergence is regulated by dynamic 
screening, the gap equation depends on po even if the frequency is small. We can also see 
that the gap scales as exp(— c/g). The collinear divergence leads to a gap equation with a 
double-log behavior. Qualitatively 
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(6) 



from which we conclude that A <~ cxp(— c/g). The approximation (6) is not sufficiently 
accurate to determine the correct value of the constant c. A more detailed analysis shows 
that the gap on the Fermi surface is given by 

A ~ 512^ 4 (2/iV / ) 5 / 2 R g- 5 cxp ("^) ■ (7) 

We should emphasize that, strictly speaking, this result contains only an estimate of the 
pre-exponent. It was recently argued that wave function renormalization and quasi-particle 
damping give 0(1) contributions to the pre-exponent which substantially reduce the gap 
[15,16]. 

For chemical potentials \i < 1 GeV, the coupling constant is not small and the appli- 
cability of perturbation theory is in doubt. If we ignore this problem and extrapolate the 
perturbative calculation to densities p ~ 5p we find gaps A ~ 100 MeV. This result may 
indeed be more reliable than the calculation on which it is based. In particular, we note that 
similar results have been obtained using realistic interactions which reproduce the chiral 
condensate at zero baryon density [9,10]. 



2.2 QCD with three flavors: Color-Flavor-Locking 

If quark matter is formed at densities several times nuclear matter density we expect the 
quark chemical potential to be larger than the strange quark mass. We therefore have to 
determine the structure of the superfluid order parameter for three quark flavors. We begin 
with the idealized situation of three degenerate flavors. From the arguments given in the 
last section we expect the order parameter to be color and flavor anti-symmetric matrix of 
the form 

4>t = (8) 
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In order to determine the precise structure of this matrix we have to extremize grand canon- 
ical potential. We find [17,18] 

A° b = A(5?5 b 3 - 5^), (9) 

which describes the color-flavor locked phase proposed in [19]. Both color and flavor sym- 
metry are completely broken. There are eight combinations of color and flavor symmetries 
that generate unbroken global symmetries. The symmetry breaking pattern is 

SU(3) L x SU(3) R xU(l)v ^ SU(3)v (10) 

This is exactly the same symmetry breaking that QCD exhibits at low density. The spec- 
trum of excitations in the color-flavor-locked (CFL) phase also looks remarkably like the 
spectrum of QCD at low density [20]. The excitations can be classified according to their 
quantum numbers under the unbroken SU(3), and by their electric charge. The modi- 
fied charge operator that generates a true symmetry of the CFL phase is given by a lin- 
ear combination of the original charge operator Q em and the color hypercharge operator 
Q = diag(— 2/3, —2/3, 1/3). Also, baryon number is only broken modulo 2/3, which 
means that one can still distinguish baryons from mesons. We find that the CFL phase 
contains an octet of Goldstone bosons associated with chiral symmetry breaking, an octet 
of vector mesons, an octet and a singlet of baryons, and a singlet Goldstone boson related 
to superfluidity. All of these states have integer charges. 

With the exception of the U(l) Goldstone boson, these states exactly match the quan- 
tum numbers of the lowest lying multiplets in QCD at low density. In addition to that, the 
presence of the U(l) Goldstone boson can also be understood. The U(l) order parameter 
is ((uds)(uds)). This order parameter has the quantum numbers of a + AA pair conden- 
sate. In Nf = 3 QCD, this is the most symmetric two nucleon channel, and a very likely 
candidate for superfluidity in nuclear matter at low to moderate density. We conclude that 
in QCD with three degenerate light flavors, there is no fundamental difference between the 
high and low density phases. This implies that a low density hyper-nuclear phase and the 
high density quark phase might be continuously connected, without an intervening phase 
transition. 



2.3 Other phases 

Color-flavor locking can be generalized to QCD with more than three flavors [17]. Chiral 
symmetry is broken for all Nf > 3, but only in the case Nf = 3 do we find the T = fi = 
pattern of chiral symmetry breaking, SU(N f ) L x SU(N f ) R — > SU(N f ) v . 

The case Nf — 1 is special [8]. In this case the order parameter is flavor-symmetric 
and the Cooper pairs carry non-zero angular momentum. For N c = 3 the spin direction 
can become aligned with the color orientation of the Cooper pair. In the color-spin locked 
phase color and rotational invariance are broken, but a diagonal 5*0(3) survives. 

If the number of colors is very large, N c — > oo, color superconductivity is suppressed 
and the ground state is very likely a chiral density wave [21]. This state is analogous to spin 
density waves in condensed matter physics. If the baryon density is large, the transition 
from color superconductivity to chiral density waves requires very large values of N c on 
the order of N c > 1000 [22]. 
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3. The Role of the Strange Quark Mass 

At baryon densities relevant to astrophysical objects distortions of the pure CFL state due 
to non-zero quark masses cannot be neglected [23-28]. This problem can be studied using 
the effective chiral theory of the CFL phase [29] (CFLxTh). This theory determines both 
the ground state and the spectrum of excitations with energies below the gap in the CFL 
phase. Using the effective theory allows us to perform systematic calculations order by 
order in the quark mass. 



3.1 CFL Chiral Theory (CFL X Th) 

For excitation energies smaller than the gap the only relevant degrees of freedom are the 
Goldstone modes associated with the breaking of chiral symmetry and baryon number. The 
interaction of the Goldstone modes is described by an effective Lagrangian of the form [29] 

2 



f 2 

C eff = ^Tr [VoSVo^t - u^E^Et] + CTr(MSt) + h.c 



(11) 

+ AiTr(MS t )Tr(MS t ) + A 2 Tr(Mtf Mtf) + A 3 Tr(MS t )Tr(M t S) + h.c] + .... 

Here £ = cxp(i(j> a A a / f n ) is the chiral field, f n is the pion decay constant and M is a 
complex mass matrix. The chiral field and the mass matrix transform as £ — > LUB) and 
M -> LMR) under chiral transformations (L,R) e SU(3) L x SU(3) R . We have sup- 
pressed the singlet fields associated with the breaking of the exact U(l)v and approximate 
U(1)a symmetries. The theory (11) looks superficially like ordinary chiral perturbation 
theory. There are, however, some important differences. Lorentz invariance is broken and 
Goldstone modes move with the velocity v n < c. The chiral expansion has the structure 

^/. 2 A 2 (|) fc (|)\sr( S t)«. (12) 

Loop graphs are suppressed by powers of d/(4irf n ). We shall see that the pion decay con- 
stant scales as f n <~ pp. As a result higher order derivative interactions are parametrically 
more important than loop diagrams with the leading order vertices. 

Further differences as compared to chiral perturbation theory in vacuum appear when the 
expansion in the quark mass is considered. The CFL phase has an approximate (Z 2 )a sym- 
metry under which M —> —M and £ — > £. This symmetry implies that the coefficients 
of mass terms that contain odd powers of M are small. The (Z 2 )a symmetry is explicitly 
broken by instantons. The coefficient C can be determined from a weak coupling instanton 
calculation and C ~ (Aq C d/pf) 8 [17,30]. BCS calculations show that the CFL phase 
undergoes a phase transition to a less symmetric phase when m 2 / (2pp) ~ A [23,24]. This 
suggests that the expansion parameter in the chiral expansion is M 2 / '(ppA). We shall see 
that this is indeed the case. However, the coefficients Ai of the quadratic terms in M turn 
out to be anomalously small 

AM 2 ~ A 2 M 2 ~ f 2 A 2 , (13) 
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compared to the naive estimate AiM 2 ~ f 2 A 2 [M 2 /(ppA)]. 

The pion decay constant f v and the coefficients can be determined using match- 
ing techniques. Matching expresses the requirement that Green functions in the effec- 
tive chiral theory and the underlying microscopic theory, QCD, agree. The pion decay 
constant is most easily determined by coupling SU(Nf)L,R gauge fields Wl,r to the 
left and right flavor currents. As usual, this amounts to replacing ordinary derivatives 
by covariant derivatives. The time component of the covariant derivative is given by 
VqS = <9q£ + iWiS — i^Wr where we have suppressed the vector index of the gauge 
fields. In the CFL vacuum S = 1 the axial gauge field Wl — Wr acquires a mass by the 
Higgs mechanism. From (1 1) we get 

C= f -f l -{W L -W R ) 2 . (14) 

The coefficients Ai^.z can be determined by computing the shift in the vacuum energy 
due to non-zero quark masses in both the chiral theory and the microscopic theory. In the 
chiral theory we have 



AS = - \A 1 (Tr(M)) 2 + A 2 Tr(M 2 ) + A 3 Tr(M)Tr(M t ) + h.c. 



(15) 



We note that different 0(M 2 ) mass terms produce distinct contributions to the vacuum 
energy. This means that the coefficients Ai can be reconstructed uniquely from the vacuum 
energy. 



3.2 High Density Effective Theory (HDET) 

In this section we shall determine the mass of the gauge field and the shift in the vacuum 
energy in the CFL phase of QCD at large baryon density. This is possible because asymp- 
totic freedom guarantees that the effective coupling is weak. The QCD Lagrangian in the 
presence of a chemical potential is given by 

£ = t(il? + ^7o) 4> ~ ^lM^jr - 4> R M^ L - \G%,G%, (16) 

where = + igA^ is the covariant derivative, M is the mass matrix and [i is the 
baryon chemical potential. If the baryon density is very large perturbative QCD calcula- 
tions can be further simplified. The main observation is that the relevant degrees of free- 
dom are particle and hole excitations in the vicinity of the Fermi surface. We shall describe 
these excitations in terms of the field ip + (i)F,x), where vf is the Fermi velocity. At tree 
level, the quark field ip can be decomposed as ip = ip + + tp_ where ip± = i(l±a-u F )^. 
To leading order in 1/pf we can eliminate the field ip_ using its equation of motion. For 
we find 

tp-,L = • Dip +tL + 7 MV>+„r) . (17) 

There is a similar equation for ip-.R. The longitudinal and transverse components of 7^ 
are defined by (70,7)11 = (jo,v(j -v)) and (7^)1. =7(u-(7^)||- To leading order in 
the lagrangian for the tp + field is given by [31-33] 
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Figure 1. This figure shows the contribution to the vacuum energy from the effective 
chemical potential terms in the high density effective theory. Fig. a) shows the diagrams 
that are matched against the (MM 1 ") 2 term. Fig. b) shows the diagram that is matched 
against the MM t EM t ME term in the chiral theory. 



C = Vl+(w • D)^ L+ - j (ipt+C^ (S ai 8 bj ~ 5 a3 5u) + h.c.) 

- T-1>l+ (W-0 2 + MM^) 4, L+ + (R^L,M~ Mt) + . . 



(18) 



with = (l,v) and i,j,... and a,b,... denote flavor and color indices. In or- 
der to perform perturbative calculations in the superconducting phase we have added a 
tree level gap term ipf- R CA ai! bjip b ^ R . In the CFL phase this term has the structure 
A aii bj = A(8 ai 5i,j — SajSbi). The magnitude of the gap A is determined order by or- 
der in perturbation theory from the requirement that the thermodynamic potential is sta- 
tionary with respect to A. With the gap term included the perturbative expansion is well 
defined. There are no additional infra-red divergences. In particular, there is no need to 
include additional gap parameters at higher order in 1/pf, such as the anti-particle gap or 
modifications of the particle gap due to non-zero quark masses [34]. 

The screening mass of the flavor gauge fields Wl.r can be determined by computing the 
corresponding polarization function in the limit qo = 0, q — > 0. We find IInn L = Ilm fl = 



-Il&f = m 2 D /A with 
equ. (14) we get [35] 



L D 



0. We find n&f 

(21 - 81og(2))p^/(367r 2 ). Matching this result against 



J IT 



21-8 



18 



log(2) (jp%_ 
\2tt 2 



(19) 



Our next task is to compute the mass dependence of the vacuum energy. To leading order 
in 1/pf there is only one operator in the high density effective theory 



£ = -- 



1 

2p~i 



(20) 



This term arises from expanding the kinetic energy of a massive fermion around p = pp. 
We note that MM t /(2p F ) and M^M/(2p F ) act as effective chemical potentials for left 
and right-handed fermions, respectively. Indeed, to leading order in the 1/pf expansion, 
the Lagrangian (18) is invariant under a time dependent flavor symmetry tpp ~ > L(t)ipL, 
ip R -> R{t)ip R where X L = MM^ j(2p F ) and X R = M t M/(2p F ) transform as left 
and right-handed flavor gauge fields. If we impose this approximate gauge symmetry on 
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the CFL chiral theory we have to include the effective chemical potentials Xl.r in the 
co variant derivative of the chiral field [27], 



/MMt\ /MtM\ 



(21) 



Xl and Xr contribute to the vacuum energy at 0(M 4 ) 
f 2 

AS = -r-^-Tr [(MM^){M^M) - (Af A/t) 2 ] . (22) 

This result can also be derived directly in the microscopic theory [27]. The corresponding 
diagrams are shown in Fig. 1. We also note that equation (22) has the expected scaling 
behavior £ ~ f 2 A 2 [M 2 /(p F A)] 2 . 

0(M 2 ) terms in the vacuum energy are generated by terms in the high density effective 
theory that are higher order in the 1/pf expansion. We recently argued that these terms 
can be determined by computing chirality violating quark-quark scattering amplitudes for 
fermions in the vicinity of the Fermi surface [34]. Feynman diagrams for (?l + <Zl — » 
qn+QR are shown in Fig. 2a. To leading order in the 1/pf expansion the chirality violating 
scattering amplitudes are independent of the scattering angle and can be represented as 
local four-fermion operators 

c = w F (^ c ^* c ^ tABCD + (^Vf ){^r R )t ACBD ) ■ 

(23) 

There are two additional terms with (L <-> R) and (Af <-» Aft). We have introduced the 
CFL eigenstates ip A defined by ipf = ip A (\ A ) m /V2, A = 0, . . . , 8. The tensors T is 
defined by 

l{Tr [\ A M(\ D f\ B M(\ c ) T ] 

- ^Tr [X A M(X D ) T ] Tr [A B M(A C ) T ] }. (24) 

The second tensor T involves both M and Aft and only contributes to terms of the form 
Tr[Af Aft] in the vacuum energy. These terms do not contain the chiral field S and there- 
fore do not contribute to the masses of Goldstone modes. We can now compute the shift in 
the vacuum energy due to the effective vertex (23). The leading contribution comes from 
the two-loop diagram shown in Fig. 2b. We find 



3A 2 

47T 2 



Tr[M]) 2 — Tr Af 2 j + (m <-» A/t) . (25) 



Using this result we can determines the coefficients A\ .2,3 in the CFL chiral theory. We 
obtain 

3A 2 

M = ~M = -j-p, A 3 = 0, (26) 
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Figure 2. Fig. a) shows the effective chirality violating four-fermion vertex in the high 
density effective theory. Fig. b) shows the corresponding contribution to the vacuum 
energy. As explained in the text, there are additional vertices which involve MM t , but 
they do not contribute to the masses of Goldstone modes. 



which agrees with the result of Son and Stephanov [35]. We also note that E ~ 
f 2 A 2 (A/p F ) [M 2 1 '(p F A)] which shows that the coefficients A4 are suppressed by 
(A/pp). The effective lagrangian (18) and (23) can also be used to compute higher order 
terms in M. The dominant (9(M 4 ) term is the effective chemical potential term equ. (22). 
Other 0(M 4 ) terms are suppressed by additional powers of (A/pp). 



3.3 Kaon Condensation 



Using the results discussed in the previous section we can compute the masses of Goldstone 
bosons in the CFL phase. In section 3.1 we argued that the expansion parameter in the 
chiral expansion of the Goldstone boson masses is 5 = m 2 /(p F A). The first term in 
this expansion comes from the 0(M 2 ) term in (11), but the coefficients A contain the 
additional small parameter e = (A/pp). In a combined expansion in S and e the O(eS) 
mass term and the 0(S 2 ) chemical potential term appear at the same order. At this order, 
the masses of the flavored Goldstone bosons are 



2 2 

mi - mi 



2 Pl 



+ 



m K± = T- 



2 2 

mi — mi 



2 Pl 



m K°.K° — T" 



mz — m j 



2pp 



+ 



4A, 

-^■{m u + m d )m s 



4A 

-j^m d (m u + m s ) 



4A . 

-j^m u (m d + m s ) 



1/2 



1/2 



1/2 



(27) 



We observe that the pion masses are not strongly affected by the effective chemical po- 
tential but the masses of the K + and K° are substantially lowered while the K~ and K° 
are pushed up. As a result the K + and K° meson become massless if m s <~ m^J^A 2 / 3 . 
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For larger values of m s the kaon modes are unstable, signaling the formation of a kaon 
condensate. 

Once kaon condensation occurs the ground state is reorganized. For simplicity, we con- 
sider the case of exact isospin symmetry m u = raj = m. Kaon condensation can be 
studied using an ansatz of the form S = cxpfmA^. The vacuum energy is 

V(a) = -fl (<^j sHa f + (m° x ) 2 (cos(a) - 1)) , (28) 

where (m {) K ) 2 = (AA/ f 2 )m Ut d(m u .d + m s ) is the 0{M 2 ) kaon mass in the limit of exact 
isospin symmetry. Minimizing the vacuum energy we obtain a = if m 2 /(2pF) < rn'x 
and cos(a) = (m^) 2 /^, 2 with [i s = m 2 /{2pp) if /j, s > mP K . We observe that the vacuum 
energy is independent of 8i,6 2} 4>- The hypercharge density is given by 

ny = ffa (l - . (29) 

We observe that within the range of validity of the effective theory, fi s < A, the hyper- 
charge density satisfies ny < Ap F /(2n 2 ). This means that the number of condensed 
kaons is bounded by the number of particles contained within a strip of width A around 
the Fermi surface. It also implies that near the unlocking transition, fi s ~ A, the CFL state 
is significantly modified. In this regime, of course, we can no longer rely on the effective 
theory and a more microscopic calculation is necessary. 

Let us summarize what we have learned about the effects of a non-zero strange quark 
mass m s . The effect of m s is controlled by the parameter m 2 /(pF A). Ifm 2 /(pi?A) <~ 1, 
color-flavor-locking breaks down and a transition to a less symmetric phase will occur. In 
the regime m 2 /(pF&) < 1 the phase structure can be established using the effective chiral 
theory of the CFL phase and dimensional analysis. We have argued that there is a new 
small scale m 2 /{pF&) ~ ( v / ?'VW^/pf) <C 1 which corresponds to the onset of kaon 
condensation. If perturbative QCD is reliable we can be more quantitative. To leading 
order in g, the critical strange quark mass for kaon condensation is 

m s | c „ t = 3.03-my 3 A 2 / 3 . (30) 

This result suggests that for values of the strange quark mass and the gap that are relevant 
to compact stars CFL matter is likely to support a kaon condensate. 



4. Conclusion: The many phases of QCD 

We would like to conclude by summarizing some of the things we have learned about the 
phase structure of QCD-like theories at finite temperature and chemical potential. We begin 
with the case of two massless flavors, Fig. 3a. If we move along the chemical potential 
axis at temperature T = 0, there is a minimum chemical potential required in order to 
introduce baryons into the system. Since nuclear matter is self -bound, this point is a first 
order transition. Along the temperature axis, the line of first order transitions eventually 
ends in a critical point: This is the endpoint of the nuclear liquid-gas phase transition. If 
we continue to increase the chemical potential, we encounter the various phases of nuclear 
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Figure 3. Schematic phase diagram of QCD at finite temperature and density. The 
figures a)-d) correspond to different numbers of massless and massive flavors, see the 
discussion in the text. 

matter at high density. Many possibilities have been discussed in the literature, and we 
have nothing to add to this discussion. At even higher chemical potential, we encounter 
the transition to quark matter and the two flavor quark superconductor. Model calculations 
suggest that this transition is first order. In the case of two massless flavors, universality 
arguments suggest, and lattice calculations support, the idea that the finite temperature zero 
chemical potential chiral phase transition is second order. In this case, the line of first order 
fi 7^ transitions would have to meet the T ^ transition at a tricritical point [36,37]. 

This tricritical point is quite remarkable, because it remains a true critical point, even 
if the quark masses are not zero, Fig. 3b. A non-zero quark mass turns the second order 
T ^ transition into a smooth crossover, but the first order /i / transition persists. 
While it is hard to predict where exactly the tricritical point is located in the phase diagram 
it may well be possible to settle the question experimentally. Heavy ion collisions at rela- 
tivistic energies produce matter under the right conditions and experimental signatures of 
the tricritical point have been suggested [38]. 

We have already discussed the phase structure of Nf = 3 QCD with massless or light 
degenerate quarks in section 2.2. We emphasized that at T = the low density, hadronic, 
phase and the high density, quark, phase might be continuously connected. On the other 
hand, there has to be a phase transition that separates the color-flavor locked phase from 
the T = n = vacuum state. This is because of the presence of a gauge invariant U(l) 
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order parameter related to superfluidity that distinguishes the two. In the case of Nf = 3 
massless flavors the finite temperature phase transition is known to be first order. We expect 
the transition from the superconducting to the normal phase at T ^ and large fi to be 
first order, too. This means that there is no tricritical point in Fig. 3c. 

The phase diagram becomes more complicated if we take into account the effects of a 
finite strange quark mass, Fig. 3d. If m 2 s /(Au) <~ A there is a phase transition between 
the Nf = 2 + 1 phase with separate pairing in the ud and ss sectors and the CFL phase 
with pairing in both ud and us as well as ds sectors. Near this phase transition we may 
encounter phases with inhomogeneous BCS (LOFF) pairing [25]. Inside the CFL phase 
kaon condensation is a possibility. 

The phase diagram shown in Fig. 3d should, at best, be considered an educated guess. 
Whether for realistic values of the quark masses there is an interlude of the the 2SC phase 
along the u ^ axis, instead of a direct transition between the CFL phase and nuclear 
matter, cannot be decided on the basis of currently available calculations. We know that 
there is at least one phase transition, because nuclear matter and the color-flavor locked 
phase are distinguished by a gauge invariant U(l) s order parameter related to strangeness. 
This, of course, is based on our belief that ordinary nuclear matter, not strange quark matter 
or hyperonic matter, is the true ground state of baryonic matter. 

Current calculations have also not conclusively answered the question whether the tran- 
sition along the T ^ axis is a smooth crossover, as indicated in figure 3d and favored by 
some lattice calculations, or whether the transition is first order, as would be the case if m s 
is sufficiently small. This is clearly an important question in connection with the existence 
of the tricritical point. 

The challenge ahead of us is to find experimental observables, both in heavy ion colli- 
sions and the observation of neutron stars, that will allow us to determine the phase diagram 
of hot and dense matter. 
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